CHAPTER 1. — 


Section 1.2 


1(a) The differential equation can be rewritten as 


dy 
5—y 


= dt. 


Integrating both sides of this equation results in — In|5 — y| = t + cı , or equivalently, 
5 —y=ce‘. Applying the initial condition y(0) = yo results in the specification of 
the constant as c = 5 — yo. Hence the solution is y(t) = 5 + (yo — 5)e™ . 


All solutions appear to converge to the equilibrium solution y(t) = 5. 


l(c). Rewrite the differential equation as 


d 
Yt. 
10 — 2y 


Integrating both sides of this equation results in — $In|10 — 2y|=t+c,or 
equivalently, 

5—y=ce ~~. Applying the initial condition y(0) = yo results in the specification of 
the constant as c = 5 — yọ. Hence the solution is y(t) = 5 + (yo —5)e~*. 


All solutions appear to converge to the equilibrium solution y(t) = 5, but at a faster rate 
than in Problem la. 


2(a). The differential equation can be rewritten as 
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d 
a A 
y—5 
Integrating both sides of this equation results in In|y — 5| = t + c , or equivalently, 
y—5= ce. Applying the initial condition y(0) = yo results in the specification of 
the constant as c = yọ — 5. Hence the solution is y(t) = 5 + (yo — 5) . 


All solutions appear to diverge from the equilibrium solution y(t) = 5. 


2(b). Rewrite the differential equation as 


dy 
= dt. 
2y—5 


Integrating both sides of this equation results in $In|2y — 5| = t + cı , or equivalently, 
2y — 5 = ce”. Applying the initial condition y(0) = yo results in the specification of 
the constant as c = 2y — 5. Hence the solution is y(t) = 2.5 + (yo — 2.5)e” . 


All solutions appear to diverge from the equilibrium solution y(t) = 2.5. 


2(c). The differential equation can be rewritten as 


dy 


= dt. 
2y — 10 


Integrating both sides of this equation results in 5ln|2y — 10| = t+ cı , or equivalently, 
y—5= ce”. Applying the initial condition y(0) = y results in the specification of 
the constant as € = yọ — 5. Hence the solution is y(t) = 5 + (yp — 5)e”. 
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All solutions appear to diverge from the equilibrium solution y(t) = 5. 


3(a). Rewrite the differential equation as 


dy 
ET = dt, 


which is valid for y 4 b /a. Integrating both sides results in —1In|b — ay| = t + cı, or 
equivalently, b — ay = ce“. Hence the general solution is y(t) = (b — ce™*)/a. 
Note that if y = b/a , then dy/dt = 0, and y(t) = b/a is an equilibrium solution. 


(b) 


(e.g. a=1 ,b=2) 


0 O02 O04 106 06 1 


(i) Asa increases, the equilibrium solution gets closer to y(t) = 0, from above. 
Furthermore, the convergence rate of all solutions, that is, a , also increases. 

(ii) As b increases, then the equilibrium solution y(t) = b/a also becomes larger. In 
this case, the convergence rate remains the same. 

(iii) If a and b both increase (but b/a = constant), then the equilibrium solution 
y(t) = b/a remains the same, but the convergence rate of all solutions increases. 


5(a). Consider the simpler equation dy, /dt = — ayı. As in the previous solutions, re- 
write the equation as 

d 

Mi _ adi. 

yı 


Integrating both sides results in y,(t) = ce™*. 


(b). Now set y(t) = y(t) + k , and substitute into the original differential equation. We 
find that 
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—ay,+0= -aly + k)+b. 


That is, — ak + b = 0, and hence k = b/a. 

(c). The general solution of the differential equation is y(t) = ce“ + b/a. This is 
exactly the form given by Eq. (17) in the text. Invoking an initial condition y(0) = yo, 
the solution may also be expressed as y(t) = b/a + (yo — b/aje™. 


6(a). The general solution is p(t) = 900 + ce", that is, p(t) = 900 + (py — 900)e"””. 
With po = 850, the specific solution becomes p(t) = 900 — 50e‘/?.. This solution is a 
decreasing exponential, and hence the time of extinction is equal to the number of 
months 

it takes, say t; , for the population to reach zero. Solving 900 — 50e" = 0, we find that 
t; = 21n(900/50) = 5.78 months. 

(b) The solution, p(t) = 900 + (po — 900)e"”, is a decreasing exponential as long as 

Po < 900. Hence 900 + (po — 900)e"”? = 0 has only one root, given by 


900 
tre = 2ln| —— |. 
i a(g) 


(c). The answer in part (b) is a general equation relating time of extinction to the value 
of 
the initial population. Setting t; = 12 months , the equation may be written as 


900 ; 
—— =e ; 
900 — pp 


which has solution p = 897.7691. Since py is the initial population, the appropriate 
answer is pọ = 898 mice. 


7(a). The general solution is p(t) = pọ e'* . Based on the discussion in the text, time t is 
measured in months. Assuming 1 month = 30 days, the hypothesis can be expressed as 
poe’ = 2p). Solving for the rate constant, r = In(2), with units of per month. 


TN/30 


(b). N days = N /30 months. The hypothesis is stated mathematically as poe = 2p 


It follows that rN /30 = In(2), and hence the rate constant is given by r = 30 In(2)/ N. 
The units are understood to be per month . 


9(a). Assuming no air resistance, with the positive direction taken as downward, 


Newton's 
Second Law can be expressed as 


in which g is the gravitational constant measured in appropriate units. The equation can 
be 
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written as dv/dt = g, with solution v(t) = gt + vo. The object is released with an 
initial 
velocity vy. 


(b). Suppose that the object is released from a height of h units above the ground. Using 
the 
fact that v = dx/dt, in which z is the downward displacement of the object, we obtain 
the 
differential equation for the displacement as dx/dt = gt + vo. With the origin placed at 
the point of release, direct integration results in x(t) = gt?/2 + vot. Based on the 
chosen 
coordinate system, the object reaches the ground when x(t) = h. Lett = T be the time 
that it takes the object to reach the ground. Then gT?/2 + vT =h. Using the 
quadratic 
formula to solve for T , 
— Ut 2gh 
T= Vo Vo + 2g , 
g 


The positive answer corresponds to the time it takes for the object to fall to the ground. 
The 

negative answer represents a previous instant at which the object could have been 
launched 

upward (with the same impact speed), only to ultimately fall downward with speed vo , 
from a height of h units above the ground. 


(c). The impact speed is calculated by substituting t = T into v(t) in part (a). That is, 


u(T) = Yw + 2gh . 


10(a,b). The general solution of the differential equation is Q(t) = ce~™. Given that 
Q(0) = 100 mg, the value of the constant is given by c = 100. Hence the amount of 
thorium-234 present at any time is given by Q(t) = 100e~™. Furthermore, based on the 
hypothesis, setting t = 1 results in 82.04 = 100e~”. Solving for the rate constant, we 
find that r = — In(82.04/100) = .19796/veek or r = .02828/day . 


(c). Let T be the time that it takes the isotope to decay to one-half of its original 
amount. 

From part (a), it follows that 50 = 100e~"", in which r = .19796/week. Taking the 
natural logarithm of both sides, we find that T = 3.5014 weeks or T = 24.51 days. 


11. The general solution of the differential equation dQ/dt = —rQ is Q(t) = Qe", 
in which Q, = Q(0) is the initial amount of the substance. Let 7 be the time that it takes 
the substance to decay to one-half of its original amount, Q,. Setting t = 7 in the 
solution, 

we have 0.5 Qo = Qye_"’. Taking the natural logarithm of both sides, it follows that 

— rr = ln(0.5) or rr = In2. 
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12. The differential equation governing the amount of radium-226 is dQ/dt = —rQ, 
with solution Q(t) = Q(0)e~—". Using the result in Problem 11, and the fact that the 
half-life r = 1620 years , the decay rate is given by r = In(2)/1620 per year. The 
amount of radium-226, after t years, is therefore Q(t) = Q(0)e~ 0042786", Let T be 
the time that it takes the isotope to decay to 3/4 of its original amount. Then setting 
t=T, 

and Q(T) = łQ(0) , we obtain 2Q(0) = Q(0)e~®-00042786T. Solving for the decay 
time, it follows that — 0.00042786 T = In(3/4) or T = 672.36 years. 


13. The solution of the differential equation, with Q(0) = 0, is 
Q(t) =CV(1—e VCR), 
As too, the exponential term vanishes, and hence the limiting value is Qz = CV. 


14(a). The accumulation rate of the chemical is (0.01)(300) grams per hour. At any 
given time t, the concentration of the chemical in the pond is Q(t) /10° grams per gallon 


Consequently, the chemical /eaves the pond at a rate of (3 x 10°*)Q(t) grams per hour. 
Hence, the rate of change of the chemical is given by 


1g 


EF 3 — 0.0003 Q(t) gm/hr. 


Since the pond is initially free of the chemical, Q(0) = 0. 


(b). The differential equation can be rewritten as 


dQ 


—__* __ = 0.0003 dt. 
10000 — Q 


Integrating both sides of the equation results in — /n|10000 — Q| = 0.0003t + C. 
Taking 

the natural logarithm of both sides gives 10000 — Q = ce~*"", Since Q(0) = 0, the 
value of the constant is c = 10000. Hence the amount of chemical in the pond at any 
time 

is Q(t) = 10000(1 — e~ °°) orams. Note that 1 year = 8760 hours. Setting 

t = 8760, the amount of chemical present after one year is Q(8760) = 9277.77 grams , 
that is, 9.27777 kilograms . 


(c). With the accumulation rate now equal to zero, the governing equation becomes 
dQ/dt = — 0.0003 Q(t) gm/hr. Resetting the time variable, we now assign the new 
initial value as Q(0) = 9277.77 grams . 


(d). The solution of the differential equation in Part (c) is Q(t) = 9277.77 e000". 


Hence, one year after the source is removed, the amount of chemical in the pond is 
Q(8760) = 670.1 grams. 
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(e). Letting t be the amount of time after the source is removed, we obtain the equation 
10 = 9277.77 e °°", Taking the natural logarithm of both sides, — 0.0003 t = 
= In(10/9277.77) or t = 22,776 hours = 2.6 years. 


(f) 


10000 
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15(a). Itis assumed that dye is no longer entering the pool. In fact, the rate at which the 
dye leaves the pool is 200 - [q (t) /60000] kg/min = 200(60/1000) (q(t) /60] gm per hour 


Hence the equation that governs the amount of dye in the pool is 


“ = —0.2q (gm/hr). 


The initial amount of dye in the pool is q(0) = 5000 grams. 


(b). The solution of the governing differential equation, with the specified initial value, 
is g(t) = 5000e-°". 


(c). The amount of dye in the pool after four hours is obtained by setting t = 4. That is, 
q(4) = 5000 e™®8 = 2246.64 grams. Since size of the pool is 60, 000 gallons , the 
concentration of the dye is 0.0374 grams/gallon . 


(d). Let T be the time that it takes to reduce the concentration level of the dye to 

0.02 grams/gallon. At that time, the amount of dye in the pool is 1, 200 grams. Using 
the answer in part (b), we have 5000 e™®™?T = 1200. Taking the natural logarithm of 
both sides of the equation results in the required time T = 7.14 hours. 


(e). Note that 0.2 = 200/1000. Consider the differential equation 
dq r 


di -1000 2" 


Here the parameter r corresponds to the flow rate, measured in gallons per minute . 
Using the same initial value, the solution is given by q(t) = 5000 e7" t0% . In order 
to determine the appropriate flow rate, set t = 4 and q = 1200. (Recall that 1200 gm of 
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—r /250 


dye has a concentration of 0.02 gm/gal). We obtain the equation 1200 = 5000 e 
Taking the natural logarithm of both sides of the equation results in the required flow rate 
r = 357 gallons per minute . 
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Section 1.3 


1. The differential equation is second order, since the highest derivative in the equation 
is of order two. The equation is linear, since the left hand side is a linear function of y 
and 

its derivatives. 


3. The differential equation is fourth order, since the highest derivative of the function y 
is of order four. The equation is also linear, since the terms containing the dependent 
variable is linear in y and its derivatives. 


4. The differential equation is first order, since the only derivative is of order one. The 
dependent variable is squared, hence the equation is nonlinear. 


5. The differential equation is second order. Furthermore, the equation is nonlinear, 
since the dependent variable y is an argument of the sine function, which is not a linear 
function. 


1g (bh) Se > yilt) = yf (t) Se" » Hence y”—y = 0. 
Also, y(t) = cosh t = yi (t) = sinht and y,/(t) = cosht. Thus y; — y, =0. 


9. y(t) = 3t + t? => y'(t) =3+ 2t. Substituting into the differential equation, we have 
t(3 + 2t) — (3t + t?) = 3t + 2t? — 3t — t? = t? . Hence the given function is a solution. 


10. y(t) = ¢/3 => yi (t) = 1/8 and yt (t) = y” (t) = yf" (t) =0. Clearly, yı (t) is 
a solution. Likewise, y,(t) = e™ + t/3 > y(t) = —e'+1/3, y(t) =e”, 

y (t) = —e™%, y?" (t) =e". Substituting into the left hand side of the equation, we 
find that e™ + 4( — e™) + 3(e™ + t/3) =e * — 4e™ + 3e™ +t =t. Hence both 


functions are solutions of the differential equation. 


11. y(t) =? > y! (t) =t-'?/2 and y” (t) = —t-*/?/4. Substituting into the left 
hand side of the equation, we have 
2° ( — 4°? /4) + 3t(¢'°/2) -tP = — 0? /2430'7/2-0'" 
=0 
Likewise, y(t) = t > yj(t) = — t” and y,’ (t) = 2t. Substituting into the left 


hand side of the differential equation, we have 2¢?(2t-*) + 3t( — t?) — t! = 4t — 
—3t'—t' =0. Hence both functions are solutions of the differential equation. 


12. y(t) =t? => y (t) = — 24° and y(t) = 6t*. Substituting into the left hand 
side of the differential equation, we have t?(6¢~*) + 5t( — 2t-°) + 4t = 6t 2 — 
—10¢?+4t? =0. Likewise, y(t) = tln t > yj(t) = t° — 2t*Int and 

ys (t) = —5t* +6t “Int. Substituting into the left hand side of the equation, we have 
t2(—5t4+ + 6t“Int) + 5t(t — 2tInt) + 4(t?2Int) = —5t? + 6t?lnt+ 
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